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ABSTRACT 

We investigate the process of synchrotron radiation from thermal electrons at semi- 
relativistic and relativistic temperatures. Wc find an analytic expression for the emis- 
sion coefficient for random magnetic fields with an accuracy significantly higher than of 
those derived previously. We also present analytic approximations to the synchrotron 
turnover frequency, treat Comptonization of self-absorbed synchrotron radiation, and 
give simple expressions for the spectral shape and the emitted power. We also consider 
modifications of the above results by bremsstrahlung. 

We then study the importance of Comptonization of thermal synchrotron radia- 
tion in compact X-ray sources. We first consider emission from hot accretion flows and 
active coronae above optically-thick accretion discs in black- hole binaries and AGNs. 
We find that for plausible values of the magnetic field strength, this radiative process is 
negligible in luminous sources, except for those with hardest X-ray spectra and stellar 
masses. Increasing the black-hole mass results in a further reduction of the maximum 
Eddington ratio due to this process. Then, X-ray spectra of intermediate-luminosity 
sources, e.g., low-luminosity AGNs, can be explained by synchrotron Comptonization 
only if they come from hot accretion flows, and X-ray spectra of very weak sources are 
always dominated by bremsstrahlung. On the other hand, synchrotron Comptoniza- 
tion can account for power-law X-ray spectra observed in the low states of sources 
around weakly-magnetized neutron stars. 

Key words: accretion, accretion discs - gamma-rays: theory - radiation mechanisms: 
thermal- X-rays: galaxies - X-rays: stars. 



1 INTRODUCTION 

The theory of cyclotron and synchrotron radiation is a well 
established part of physics. However, there still remain un- 
certainties about the accuracy and the range of applicabil- 
ity of some analytic formulae describing the emission. One 
important example of such uncertainties concerns the spec- 
tra of synchrotron emission from mildly-relativistic and rel- 
ativistic thermal plasma, in which case numerous studies 
devoted to this field, e.g., Jones & Hardee (1979), Petrosian 
(1981, hereafter P81), Takahara & Tsuruta (1982), Petrosian 
& McTiernan (1983) , Robinson & Melrose (1984), Mahade- 
van, Narayan & Yi (1996, hereafter MNY96), yielded results 
not entirely consistent with each other. 

Precise determination of spectra from the thermal syn- 
chrotron process is of key significance for studies of emission 
from accretion flows onto black holes and neutron stars. 
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Although direct, optically-thin, thermal synchrotron emis- 
sion is rarely observable, accurate optically-thin spectra are 
necessary to determine the turnover frequency, below which 
the plasma becomes optically-thick. The resulting partially- 
absorbed spectrum is then observable in some cases. Fur- 
thermore, photons from that spectrum provide a supply of 
seed photons for Comptonization in the plasma, which pro- 
cess gives rise to observable power-law spectra with high- 
energy cutoffs. On the other hand, blackbody radiation emit- 
ted by optically-thick accretion discs and other forms of 
optically-thick matter may constitute a competing supply 
of seed photons. In addition, bremsstrahlung radiation is 
always emitted by a hot plasma. 

In this work, we first consider formulae for the syn- 
chrotron emission coefficient of a thermal plasma (Section 
We clarify the accuracy and the range of applicability 
of previously derived formulae as well as propose our own 
expressions. We concentrate on mildly relativistic and rel- 
ativistic plasmas, as exhaustive studies of non-relativistic 
thermal cyclotron emission exist (Chanmugam et al. 1989 
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and references therein). We then calculate the turnover fre- 
quency and consider effects of bremsstrahlung emission and 
self-absorption. Second, we investigate Comptonization of 
the synchrotron radiation (hereafter abbreviated as the CS 
process) and present convenient formulae for the resulting 
spectra and luminosities (Section ^ . In Section we apply 
our results to two main geometries of accretion flows onto a 
black hole: a hot, two-temperature, optically thin disc, and 
active regions above a cold disc (i.e., a patchy corona). Those 
two models, for various plausible prescriptions for the mag- 
netic field strength, are then compared with spectral data 
from black-hole binaries and AGNs. Finally, in Section^ we 
study the importance of the CS proces for formation of spec- 
tra of accreting neutron stars with weak/moderate magnetic 
fields. 



2 SYNCHROTRON RADIATION FROM 
THERMAL PLASMAS 

2.1 Emission of a single electron 

Let us consider an electron (with charge e) moving in a 
uniform magnetic field, B, at a velocity, (3 = v/c. Let ^ be 
the angle between v and B, and ■& be the angle between 
B and the direction towards the observer. Then, the cyclo- 
synchrotron power per unit frequency and unit solid angle 
in the observer frame and in cgs units is given by (e.g. Bekefi 
1966; Pacholczyk 1970) 



(1) 



n = l 



'COSl9 - /3C0SC ^ 2, , , r,2 ■ 2tT'2/ ^ 

' J„(z) +13 sm [z) 



where 

nvc , „ ^ „s 1/7/3 sin I? sin f ,„, 

yn = — - - u(l - f3cos£cosi)), z = — ^, (2) 

7 

Uc = eB /2nmcC is the cyclotron frequency, 7 = (1 — 0^)~^^'^ 
is the Lorentz factor, J„ is a Bessel function of order n, and 
me is the electron mass. 

Note that an additional factor of (1 — /3/icosi9)~^ due 
to the Doppler effect should have appeared in the formal 
derivation of equation (0). However, this term disappears 
in the case of an electron moving chaotically, as, e.g., in a 
thermal plasma. For detailed discussion, see Scheuer (1968), 
Rybicki & Lightman (1979, section 6.7), Pacholczyk (1970, 
section 3) and references in these works. 



2.2 Synchrotron emission coefficient in a thermal 
plasma 

The emission coefficient, jv{'&), of a thermal plasma at a 
temperature, T, and with a uniform magnetic field, B, can 
be obtained by integrating the rate of equation (|l|) over a 
relativistic Maxwellian electron distribution, 



nc 7(7^ - 1) 



e 



1/2 



■ exp 



(3) 



A'2(i/e) 

where O = kT/rricC? is the dimensionless plasma tempera- 
ture, and Tic is the electron density. K2 is a modified Bessel 



function, which can be approximated by 




1 J- ise I lose^ 



, ln(2e) + 3/4-7E 

+ se^ 



+ 



ln(2e)+0.95 



-1/e 



(4) 

e ^ 0.65; 
e ^ 0.65, 



where 7e ~ 0.5772 is Euler's constant and the last coeffi- 
cients in the series have been adjusted to achieve a relative 
error ^ 0.0008. 

The integral for the emission coefficient is then, 

<il^M"/) J d^ir?^(i?,^,7), (5) 

where fi — cos 5. This integration is relatively difficult to 
carry out due to the complicated form of the integrand and 
the presence of J„ in rji,. In particular, standard numerical 
methods of computing J„ become very inefficient for n 3> 
1. P81 obtained an approximate solution by replacing the 
summation over n by integration over v in equation (|l|) and 
then using the principal term of an asymptotic expansion 
of Jn{z), where ^ 2: ^ n/3 < n [see equation (9.3.7) of 
Abramowitz & Stegun (1970, hereafter AS70)]. The region of 
validity of that approximation of J„ is given by 7^ <^ I'/vc- 
The integration over ^ and 7 by the method of the steepest 
descent with some additional approximations yields 



30^-2(1/6) 



exp 



1/3' 



(6) 



where v = u/uc&^. We see that the emission coefficient be- 
comes very small both at large frequencies, and at viewing 
angles, t?, significantly different from n/2. 

We note that equation (26) in P81, which corresponds 
to the formula (^), is too small by a factor of 2 due to the 
adopted normalization of the electron distribution [equation 
(23) in P81] being also twice too small. On the other hand, 
Takahara & Tsuruta (1982) obtain a formula [equation (2.9) 
in their paper] for the absorption coefficient for the per- 
pendicular polarization of synchrotron radiation, , which 
agrees with our expression (^) for the emission coefficient at 
i9 — n/2. (They also compute the coefficient for the parallel 
polarization, a|l, which, however, is negligibly small com- 
pared to the dominant coefficient of the perpendicular po- 
larization.) We note that the relation between the absorption 
and emission coefficients is given by Kirchhoff's law, which, 
in the case of polarized radiation, contains the source func- 
tion for each polarization separately, i.e., bJ] = B^^r = 3^/2 
(see, e.g. Chanmugam et al. 1989). The coefficients for polar- 
ized radiation are then related to the coefficients including 
both polarizations by = jjl -I- j^, Oi, = (all -I- a'^)/2. We 
also note that the emission coefficient computed by Jones & 
Hardee (1979) for 2> 1 is larger than the corresponding 
limit of equation (|^ by a factor of 2^^'^, which discrepancy 
we have not been able to explain. 

In our numerical calculations, we used an approxima- 
tion of Jn{z) in terms of the Airy function [equation (9.3.6) 
in AS70], which maximum relative error is ~ 0.08 at n = 1, 
z = 0, but it rapidly decreases for larger values of n, z. If the 
argument of the Airy function is < 2.25, we approximate it 
by its power series [equation (10.4.3) in AS70] up to the 13th 
power. Otherwise, we use the approximation to J„ of P81 
described above, in which case the maximum error is ~ 0.02. 
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Figure 1. (a) Contour plots of the ratio of tlie approximated 
emission coefficient of equation (^) to the exact numerical value 
(at ■& = Tc/2). (b) The same for equations 



Typical relative accuracy of the resulting approximation to 
Jn is then < 0.01. The accuracy can be further increased by 
adding the first-order correction to the principal term of the 
asymptotic expansion of J„, see equation (9.3.7) in AS70. 
The derivative, J^, is calculated with the second expression 
in equation (9.1.30) of AS70. Note that the above method 
is more accurate, but also more complicated, than a related 
method given by Wind & HiU (1971). 

We have tested the accuracy of formula ^ compared 
to the results of numerical integration (which themselves are 
in a very good agreement with numerical results tabulated 
for ^ 0.1 by Chanmugam et al. 1989) of equation ^ at 
= 7t/2 for 10 < v/fc 5? 1000 and 10"^ sS 9 10. These 
ranges include parameters most relevant for compact sources 
(e.g. Takahara & Tsuruta 1982; Zdziarski 1986, hereafter 
Z86; Narayan & Yi 1995, hereafter NY95), and we hereafter 
use them while testing other synchrotron rates. Figure nk 
shows the results for O > 0.02. We see that equation M 
overestimates the actual value of by orders of magnitude 
at low values of B, and should not be used at all for O 
0.1. The relative accuracy strongly depends on Q (being 
best at © ~ 0.4) and it improves with increasing v/vc- In 
the relativistic limit, O 2> 1, the accuracy of formula (^ 
depends on v only. 

We thus see that the accuracy of the emission coefficient 
(P) is not satisfactory for detailed modelling of astrophysi- 
cal plasmas, as also found by MNY96. We suggest using a 
significantly more accurate expression [eq. (11) of P81 with 



additional approximations of eq. (25) of P81 and eq. (20) of 
Petrosian & McTiernan (1983)], 
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Here 
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and 



70 = < 



1+ 



2 i/c 



e 0.08; 
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e ^ 0.08, 



where 70 is the Lorentz factor of the saddle point of the inte- 
gral (|^ over 7 (i.e. of those thermal electrons that contribute 
most to the emission at v), corresponding to the minimum 
of 720(7)2(1?, 7)77. Note that approximations (^ and (|lo|) 
have a small discontinuity at ~ 0.08. 

Figure |^b shows the relative accuracy of the approxi- 
mation of equations (^)-(|l^. We see it is much more ac- 
curate than expression (n) throughout the test range, with 
the minimum relative error at 0.01 ^ 0.2. Also, the 
relative accuracy of this approximation varies with Q and 
v/uc much slower than that of equation (^. At ~ 10~^ 
the approximation (0)-(|l^) matches well the nonrelativistic 
emission coefficient given by equations (14)-(15) of Trub- 
nikov (1958). 

Our expression (0)-(0) can be compared with that of 
Robinson & Melrose (1984), who have also used a method 
based on that of P81 with accuracy improved with respect to 
that paper. We have found that their expression, in the form 
of Dulk (1985) and with the non-relativistic Maxwellian 
replaced by the relativistic one, is slightly more accurate, 
but also more complicated, than expression ([^-(|lo|). We 
also note that a relativistic generalization of the results 
of Robinson & Melrose (1984) by Hartmann, Woosley & 
Arons (1988) appears incorrect. Namely, their expression 



(Bl) should be multiplied by a factor 271(79 
their notation). 



-1/2 



(in 



2.3 Angle-averaged emission coefficient 

The emission coefficients derived above are appropriate for 
a plasma in a uniform magnetic field. However, in typical as- 
trophysical conditions, we expect either to deal with emis- 
sion from regions where magnetic fields are chaotic or to 
observe radiation being a sum of contributions from a num- 
ber of regions with different orientations of magnetic field. 
Therefore, we would like to find the emission coefficient av- 
eraged over magnetic field directions (or, equivalently, over 
the viewing angle), = ei^/47t, where is the emitted 
power per unit volume and frequency, and 
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>(i9)smi? dj?. 



(11) 



To obtain an expression for j^, we integrate equation 
(^) over 1? using the method of the steepest descent with the 
saddle point ai ^ — n/2. Here, we treat Z as the fast varying 
(with i9) part of the integrand The resulting expression 



2ic 



5' ln2(t?,7o) 



(12) 



to be evaluated at = n/2. 

On the other hand, the asymptotic emission coefficient 
(h) can be integrated over as done by MNY96. This yields. 



35/6cA'2(l/e)t;i/6 



a(0, v) exp 



(t) 



1/3' 



(13) 



where the correction factor, a, represents the ratio of the 
exact emission coefficient to that obtained by integration of 
equation (^. This factor can be calculated from the ratio of 
equations (33) to (31) in MNY96 with fitting coefficients of 
their Table 1. 

We have tested the accuracy of the formulae ( |l^ ) and 
(p^, the latter both with and without the tabulated correc- 
tions, again in the range 10""^ ^5 ^ 10 and 10 ^u/uc ^ 
1000. Figure ^a shows the relative error of equation (^^. We 
see that the relative error is typically ^ 30 per cent except 
for e > 1 with the best accuracy at 0.02 <, Q <, 0.4. 

The expression (^|) with the tabulated corrections of 
MNY96 matches very well the numerical results (the dif- 
ference is less than a few per cent), with the exception of 
the case of = 0.084, where the corrections in Table 1 of 
MNY96 appear misprinted as the obtained values are a few 
times too small and the error increases with frequency. 

The accuracy of formula ( ^) with a = 1 is shown in Fig- 
ure We see that it provides an estimate correct within 
a factor of ~ 3 or better only for O > 0.1 and the best 
accuracy is obtained at O 0.6 This formula strongly over- 
estimates the correct result for O 0.1. At relativistic tem- 
peratures, 3> 1, the relative accuracy of equation ( |l3| ) 
with a = 1 depends on v only. It underestimates somewhat 
the actual emission coefficient, and its integration over v 
yields a value of the total emitted power too low by a factor 
of 3^/^357t/2* ~ 0.744 with respect to the actual power, 



B 

e = 166 ncfJTC — , O S> 1, 

87t 



(14) 



where ctt is the Thomson cross section. Then, equation ( p^ 
with a = 1 and divided by 0.744 represents an approxima- 
tion to the ultra-relativistic thermal synchrotron emission 
coefficient satisfying the total-power constraint (|l^). 



2.4 Synchrotron self-absorption 

The synchrotron radiation is self-absorbed by electrons up to 
the turnover frequency, Ut, above which the plasma becomes 
optically thin to the synchrotron radiation, i.e.. 



R = l, 



(15) 



where R is the characteristic size of the plasma. Below ft, 
the observed spectrum has the blackbody form, typically in 
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Figure 2. (a) Contour plots of the ratio of the approximated 
emission coefficient j|l2|) to the correct numerical value, (b) The 
same for equation (^IsK^with a = 1. 



the Rayleigh- Jeans limit. In this limit and averaging over 
angles, Kirchhoff 's law implies. 



-R 



(16) 



This can be solved numerically. On the other hand, Zdziarski 
et al. (1998) give the solution using equation (|l^). 



_ 343 C 
Q^uc 36 In " ' 



_3_ 

79 



7TTTCiexp(l/0) 



.(17) 



where Xc = hvc/nieC is a dimensionless cyclotron frequency, 
tt = UcO-tR is the Thomson optical depth of the plasma, 
and Of is the fine-structure constant. Note that i^t depends 
on rio and R only through their product, or equivalently, 
Tt. Typically, the correction factor, a, is a slowly- varying 
function of frequency (MNY96) and then equation ( p^ ) for 
ft is explicit. We also see that the dependence of Vt on a 
is only logarithmic, and thus the accuracy of determining 
the synchrotron emission coefhcient only weakly affects the 
turnover frequency. 

Mahadevan (1997) used a slightly different method of 
determining ft- Namely, he equated the total flux of the 
Rayleigh- Jeans emission to that in the optically-thin syn- 
chrotron one. This corresponds to setting the right-hand side 
of equation jl^ ) to 3/4 instead of 1, which changes slightly 
the resulting value of Ut as compared to equation (|l7|). 
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Figure 3. The ratio of the turnover frequency calculated from 
equation with a = 1 to the accurate numerical result as a 
function of temperature. 



We have compared the approximated formula ( |l7| ) with 
a = 1 with results of numerical calculations. It turns out 
that the error of equation is almost independent of the 
values of the parameters other then O. Figure |^ compares 
the results of the two methods in the temperature range 
0.03 ^0^1 for a source with tt ~ 1, R = 10^ cm and B = 
10^ G. (In this parameter range, the effect of bremsstrahlung 
can be neglected, see below.) We see that for 0.1 6 ^ 
1 the discrepancy does not exceed 20 per cent but grows 
rapidly for lower temperatures. We have also found a power- 



law approximation for Xt 



xt ~ 2.6 X 10"'' — 



e 

0.2 



hut/nieC 
B 



0.05 



10«G 



(18) 



which is accurate to 30 per cent for 0.05 ^ ^ 0.4, 
0.3 ^ TT ^ 3, 10 G ^ B ^ 10* G. Typical values of xt/x^ at 
mildly-relativistic temperatures are in the range of 10-10"^. 



2.5 Effect of bremsstrahlung on the turnover 
frequency 

We note that the above formulae for the turnover fre- 
quency do not take into account bremsstrahlung emission 
and self-absorption. However, for some combinations of 0, 
Tie and B, bremsstrahlung emission can be comparable to 
or stronger than the synchrotron one and then the turnover 
frequency determined neglecting bremsstrahlung will be in- 
correct. Then, the absorption coefficient, a^^ , should include 
a contribution from bremsstrahlung, . Note that since 
oc n^R, Vt is no more solely a function of tt whenever 
bremsstrahlung is important. 

The effect of bremsstrahlung on the turnover frequency 
can be checked by computing Vt without taking into ac- 
count bremsstrahlung, and then computing the emission co- 
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Figure 4. Example dependences of the turnover frequency on 
(a) temperature and (b) magnetic field strength. The dot-dashed 
and dashed curves give vt calculated with the analytic approx- 
imation of equation ([l7[), and numerically using equation 
respectively, with bremsstrahlung neglected in both cases. The 
solid curves give the numerical results with bremsstrahlung taken 
into account. 



efficient, , includi ng b oth processes. As long as ji,^ 3> j^^ , 



2.4 can be used. We hereafter use formu- 



equations in Section 
lae for the free- free emission coefficient of Svensson (1984). 

We have calculated the turnover frequency as a function 
of and B for the remaining parameters fixed, see Figures 

b, respectively. As expected, Ut is determined mostly 
by bremsstrahlung at low temperatures and weak magnetic 
fields. However, though bremsstrahlung can have a negligi- 
ble effect on the value of Vt , the total emitted power can still 
be dominated by it, see Section ^ below, and, e.g. Figure 3 
in Z86. On the other hand, if bremsstrahlung self-absorption 
dominates over the synchrotron one, bremsstrahlung emis- 
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sion will also dominate over the CS emission, except in the 
case of spectra of the latter emission being harder than the 
bremsstrahlung spectrum. 



3 COMPTONIZATION OF SYNCHROTRON 
PHOTONS 

Synchrotron photons produced in a plasma will, in general, 
be Compton scattered by the plasma electrons. In the case 
of a thermal plasma, the synchrotron spectrum is usually 
self-absorbed up to a high value of v/uc (Section 2.4). In 
that case, the self-absorbed synchrotron spectrum is rather 
narrow, consisting of a hard Rayleigh- Jeans spectrum at 
V ^ (in which regime scattering can be neglected, r 3> tt, 
T 1), and a fast-declining tail of the optically-thin syn- 
chrotron spectrum &t v ^ Vt ^ Vc (where absorption can be 
neglected, r <g tt, t ^ 1). Photons in that spectrum then 
serve as seed photons for Comptonization. For that process, 
it is usually sufficient to treat the self-absorbed synchrotron 
photons as monochromatic at Vt (see, e.g., Monte Carlo sim- 
ulations in Z86). Also, for magnetic fields, B ^ 10® G, and 
electron temperatures, kT > 10 keV, expected in compact 
sources, hvt <C kT, and photons &t v ^ i^i gain energy in 
the scattering process. 

To treat thermal Comptonization, we use the method of 
Zdziarski (1985, hereafter Z85), as applied to thermal syn- 
chrotron sources in Z86. This method gives an approximate 
form of Comptonized spectra, reproducing relatively accu- 
rately the energy spectral index, a, and the total power in 
the scattered spectrum (see comparison with Monte Carlo 
results in Z86). On the other hand, it gives a relatively inac- 
curate shape of the high-energy cutoff of the scattered spec- 
trum (see, e.g., Poutanen & Svensson 1996), which, however, 
is of negligible importance for our applications. 

We consider a homogeneous and isotropic source char- 
acterized by tt, Q and B. The fiux in scattered photons can 
be approximated as a sum of a cut-off power law and a Wien 
component, 



i^c(x)=iVp(|) "e-^/®-fiVw 



(5)' 



(19) 



where x = hv/mc(? is a dimensionless photon energy, and 
A'p, A'w are normalizations of the power-law and Wien com- 
ponents, respectively. 

For Tt 2, a can be approximately given by, 

a = -^, A^l + (^)^l+4e + 16e^ (20) 

where A is the average photon energy amplification per scat- 
tering and Pbc is the scattering probability averaged over the 
source volume. Note that a > in general. In spherical ge- 
ometry (Osterbrock 1974), 



Psc = 1 



8r| L 

4 ' „ 



1 



4t'p ' 



-l + e 
TT < 1; 

TT > 1, 



2tt 



(2tt + 1)] 



(21) 



where tt is the Thomson optical depth along the radius. We 
have found a power-law approximation to a(rT, 6) valid at 
0.5 TT 2 and e 0.4, 



(22) 



See Zdziarski et al. (1994) for Psc in the slab geometry. In 
general, different prescriptions for a(r, Q) can be used, e.g., 
in order to account for a different geometry or to increase 
the accuracy, without affecting the validity of the results 
presented below in this section. On the other hand, some 
equations in Section ^ utilize equation (^2|), and using a 
different a(T, 0) would somewhat change those formulae. 

Note that at tt <C 1, scattering profiles corresponding 
to subsequent scatterings are visible in the Compton spec- 
trum, and a power law is no longer a good approximation to 
the shape of the spectrum below the cutoff (e.g.. Fig. 2b in 
Z86). Still, the power- law description presented here gives 
the shape averaged over the scattering orders, and an inte- 
gral over that approximate form gives a fair approximation 
to the total luminosity (see, e.g., Monte Carlo simulations in 
Z86). On the other hand, at tt ^ 2-3, Comptonization can 
be described by means of a kinetic, Fokker-Planck, equation 
(Sunyaev & Titarchuk 1980; see e.g. Lightman & Zdziarski 
1987 for relativistic and Iow-tt corrections). The two solu- 
tions can be matched at tt ^ 2 (Z85). The remaining results 
presented in this Section can still be used at tt 2-3, ex- 
cept for a different prescription for a. In particular, the ratio 
of Nw/Np given by Z85, 



(23) 



A^w _ r(Q) 

Np r(2a + 3) 

where F is Euler's gamma function, and which uses a result 
of Sunyev & Titarchuk (1980), is a good approximation in 
both the optically-thin and optically-thick regimes (see, e.g., 
comparison with Monte Carlo results in Z86). 

We need then to normalize the fiux in the Compton- 
scattered spectrum, Fc{x), with respect to that in the 
Rayleigh- Jeans spectrum, -Frj(x). In general, the self- 
absorbed synchrotron spectrum at its peak around xt (= 
hvt / mc(?) is somewhat above an extrapolation of the power- 
law component of the Compton spectrum, Fc{x) (which fol- 
lows from photon conservation, Z85) . Z86 finds the following 
phenomenological relation providing a good approximation 
to the relative normalization. 



1 + {2Qf 
1 + 10(29)2 



(24) 



This is valid for a;t <C O, which condition we assume here- 
after. The flux of the Rayleigh- Jeans spectrum is given by, 



Fbj{x) = 7T/Rj(a;), Irj{x) = ^"^"3 ^ a:^. 



(25) 



where Jrj is the specific intensity and Ac — 2.426 x 10 
cm is the electron Compton wavelength. Then, 

Np ^niplRj{xt)ixt/er. (26) 

The luminosity due to the CS emission is then given by. 



Lcs = A 



dx Frj (x) + 



dx Fc (x) 



(27) 



where A is the source area, and which can be integrated to. 



Lcs = A 



27tmcC^6a;? 



1 + 3V(| 
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6r(Q)Pec 



r(2Q + 3) 



(28) 



where the incomplete gamma function is well approximated 
for a;t <C 6 by, 



'('-"•I) 



In 



e 



- 7e, 
1 /'e^"-l 



r(i-a) 



+ 



|a- 1[ < 10" 
|a - 2| < 10- 
a > 2.99; 
otherwise. 



(29) 



At xtlQ — 0.01, the maximum relative error of this ap- 
proximation of ~ 0.02 occurs around a ~ 3. The relative 
error declines rapidly at lower values of Xt/O and a; e.g., 
it is < 0.001 at Xt/Q = 0.01 and a sC 2.9, and < 0.002 at 
Xt/Q = 10""^ at any a. 

In an astrophysically important case of a ~ 0.4-0.9 (e.g. 
Gierliiiski et al. 1997; Zdziarski, Lubiriski & Smith 1999), the 
main contribution to the luminosity comes from the high- 
energy end of the spectrum, and the Wien component is 
relatively unimportant. Then, we obtain an approximation 
valid within a factor of 2 at Xt/Q 10~^, 



Lcs 



(1 - a)Xl 



(30) 



where we assumed a spherical geometry. On the other hand, 
for soft spectra, with a ^ 1.1, we get an approximation valid 
within ~ 30 per cent, 

Lcs ~ 7^ 1 + — ^ a^t e. (31) 



3\l 



For order-of-magnitude estimates, we can then substitute xt 
of equation (Flq) in equations (pCl)-(pl|). This yields Lcs oc 

D2o3.78-0.11<t^l.82+0.91a ^U^tk^rCiG ^2.73 r^^ n A <^ „. 



for OA<,a<, 



0.9 and a ^ 1.1, respectively. 

We stress that the above formulae for the Comptoniza- 
tion spectral shape and the corresponding expressions for 
the luminosity are mutually consistent. This is much prefer- 
able to computing Comptonization luminosity independently 
of its corresponding spectrum by multiplying the power in 
a self-absorbed synchrotron spectrum by an amplification 
factor of Comptonization (e.g. of Dermer, Liang & Canfield 
1991), as often done in studies of advective flows. 



4 APPLICATIONS TO ACCRETING BLACK 
HOLES 

Figure ^ shows an example of the CS spectrum for parame- 
ters typical to low-luminosity AGNs. The spectrum has the 
Rayleigh- Jeans shape below the turnover frequency, it is a 
power law due to Comptonization above that frequency, and 
then it has a thermal high-energy cutoff. We also show the 
contribution due to Comptonized bremsstrahlung, which is 
moderately important for the chosen parameters. We have 
computed the latter spectrum by treating the spectrum of 
equation ( |l9| ) as Green's function for Comptonization and 
then integrating over the seed spectrum of optically-thin 
bremsstrahlung. The spectrum is computed assuming spher- 
ical geometry. More examples of spectra from magnetized 
plasmas are given, e.g. in Z86. 
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Figure 5. An example of the spectrum of a magnetized plasma, 
for _R = 9 X 10^3 cm, S = 2.2 X 10^ G, tt = 0.47 and = 0.4. 
The dotted, short-dashed and long-dashed curves represent spec- 
tra due to synchrotron emission, its Comptonization and Comp- 
tonized bremsstrahlung, respectively. The solid curve represents 
the sum of those components. This spectrum also represents a 
model of a hot accretion flow fitted to the X-ray data for NGC 
4258 (bow-tie, Makishima et al. 1994), see Section 4.1. 



For given O, the Rayleigh- Jeans part of the spec- 
trum is independent of B, and its spectral luminosity is 
L^i(y) oc _R^. On the other hand, the CS spectrum obeys 
Lcsiy) oc E?v^^°' [equation (p6[)], i.e., its normalization in- 
creases quickly with increasing B via the dependence of the 
turnover frequency on B. When self-absorption is dominated 
by the synchrotron process, Vt is roughly oc B with no ex- 
plicit dependence on _R, see equation ([18|). The total luminos- 
ity, Z/cs, follows the same dependence if a 1, see equation 
(po|), and Lcs oc B?ul for q <; 1, see equation (^). The 
shape of the Comptonized bremsstrahlung spectrum is al- 
most independent of B (except for a very weak dependence 
of Comptonization on Vt). 

The bremsstrahlung luminosity can be constrained from 
below as 



LcB 



> 1.2 X io"'^e^/V|r- fa 1.2 X io"*Q"^/^e"V, 



(32) 



where r = R/Rg, Rg = GM/(? is the gravitational ra- 
dius, Le = 4:niicGMmpc/aT ~ 1.5 x 10"^*m erg s"""^ is the 
Eddington luminosity, /ic — 2/(1-1- X) is the mean elec- 
tron molecular weight, and X (~ 0.7) is the H mass frac- 
tion, and m = M/Mq. The right-hand-side expressions ne- 
glect Comptonization and relativistic corrections, which ef- 
fects will increase the actual Lcb- The last expression uses 
the approximation of equation ( p^ ) to a(rT,0). Note that 
Lcb represents a strict lower limit to the luminosity of a 
source with given tt, O and size. Thus, the luminosity of 
weak sources is dominated by Lcb, with Lcs being negli- 
gible. Note, however, that usually the plasma parameters 
depend on luminosity, which dependence should be taken 
into account when determining the luminosity below which 
bremsstrahlung dominates. 

Figure ^ shows some example dependences of the total 
L on B at e = 0.05-0.4 for i? = 4 x lO'^ cm and 4 x lO" cm 
(corresponding to ~ 30-Rg for IOM0 and lO^M© black-hole 
mass, respectively). The flat parts at low values of B cor- 
respond to the dominant bremsstrahlung, with its relative 
importance increasing with decreasing O. 

Hereafter, we use the symbols L and rj for the total 
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Figure 6. Example dependences of the plasma luminosity on 
B for (a) _R = 4 X lO'^ cm and (b) R = 4 X 10^"' cm for three 
values of temperature. The flat parts correspond to dominant 
bremsstrahlung, and dashed curves give the luminosity in the CS 
component alone. 



luminosity of a hot plasma in a source, including all ra- 
diative processes, and for the corresponding Eddington ra- 
tio, L/Le, respectively. This L then corresponds to the one 
observed (excluding components not originating in the hot 
plasma, e.g., a disc blackbody). We then compare observed 
values of L with the CS luminosity assuming some specific 
prescriptions for the magnetic field strength in an accre- 
tion fiow. Typically, the maximum possible magnetic field 
strength corresponds to equipartition between the pressure 
or energy density of the field and of the gas and radiation 
in the source (Galeev, Rosner & Vaiana 1979). 



4.1 Two-temperature accretion flows 

Here, we consider the CS emission from optically thin, two- 
temperature accretion fiow. In inner parts of hot accretion 
flows, the ion temperature, Ti, is typically much higher than 
the electron temperature (e.g. Shapiro, Lightman & Eardley 
1976; NY95), and the ion energy density is much higher than 
that of both electrons and radiation. Then, energy density 
equipartition corresponds to 



87t 2 /ii ' 



(33) 



where = 4/(1-1- 3X) is the mean ion molecular weight. 
We note that pressure equipartition would result is a slightly 
different condition, and that the magnetic-field pressure, 
B'^/8n, is sometimes assumed to equal B^/24tc (e.g., Ma- 
hadevan 1997). 

Maximum possible ion temperatures are sub-virial, and 
detailed accretion flow models show kTi ~ SnipC? /r, where 
5 <C 1 (Shapiro et al. 1976; NY95). In particular, 5 constant 
with r and dependent on the flow parameters is obtained 
in the self-similar advection-dominated solution of NY95. 
Close to the maximum possible accretion rate of the hot 
flow, 5 ~ 0.2 is a typical value [see equation (2.16) in NY95], 
which we adopt in numerical calculations below. However, 
the self-similar solution breaks down in the inner flow, which 
results in values of Ti much lower than that given by the self- 
similar solution, e.g., by a factor of ~ 10 at r = 6 (Chen, 
Abramowicz & Lasota 1997). Thus, we can obtain an up- 
per limit on the equipartition magnetic field strength by 
assuming that it follows the self-similar solution above some 
radius, ro ~ 20, while it remains constant at r < ro, 

1/2 



B ~ 7.3 X 10" G 



(^) 



lAo, 

l/r. 



r < ro; 
r > ro. 



(34) 



This typically gives B ^ 10^ G and lO"* G m mner regions 
of binary sources and AGNs, respectively. 

Then, the CS luminosity can be constrained from above 
by a sum of the contribution from the inner region (assumed 
here to be spherical) with radius ro, and a contribution from 
the outer region r > ro, obtained here by radial integration 
all ro to infinity (for which we assumed a planar geome- 
try). We use here the upper limit on B of equation (M) and 
assume a and 6 constant through the fiow. The latter as- 
sumption maximizes Lcs of given a and since both tt 
and <d will, in fact, decrease with r in an accretion fiow. 
(We also note that the maximum of dissipation per Inr 
in the Schwarzschild metric occurs at r ~ ro ~ 20.) For 
0.4 Q ^ 0.9, the resulting upper limit is 

Lcs _ 43 X 10-2.83^g0.91+0.46^g)2.64-0.68^^ 
Ti" ~ (mr2)0-46o^-0.09(l _ „)Q,0.63a + 1.26 ^ 

f 1 for inner region; 

\ (0.9lQf — 0.18)^^ for outer region, 

where we used equations (|^, ^ and ~ 0.5. The 

total Lcs corresponds to the sum of the contributions from 
the two regions. In the case of a ^ 1.1, we have an upper 
limit on the luminosity of 

Lcs . 0.0215i-^^ei-»'= 3(/p 



(35) 



^0.37^0.73 



Qi.oo Y a — 1, 
1 for inner region; 
1.37 for outer region, 



(36) 
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which is based on equation (|3l| ). 

We see from equations (|3^-(^6|) that the predicted Ed- 
dington ratio decreases with the mass. Then, if the CS pro- 
cess dominates, the predicted X-ray spectra would harden 
with increasing M at a given Lcs/Le- In fact, the opposite 
trend is observed; black-hole binaries have X-ray spectra 
harder on average than those of AGNs (e.g. Zdziarski et al. 
1999), implying that either the CS process does not domi- 
nate the X-ray spectra of those objects or (less likely) the 
Eddington ratio is much lower on average in AGNs than in 
black-hole binaries. 

We note that this conclusion differs from a statement 
in NY95 that hot accretion flows with CS cooling are effec- 
tively mass-invariant. In particular, we have been unable to 
explain their Fig. 4, in which flows onto black holes have ra- 
dial dependences of both the temperature at the critical (i.e., 
maximum possible) accretion rate and that rate itself (in Ed- 
dington units) virtually identical for m — 10 and m — 10*. 
This would then imply very similar X-ray spectra in both 
cases. On the other hand, we have found a good agreement 
of our results with the dependences on m in the results of 
Mahadevan (1997). 

We find that, under conditions typical to compact 
sources, the Comptonized bremsstrahlung luminosity, Lcb, 
is often comparable to or higher than Lcs- Figure M shows 
the results of numerical calculations of the Eddington ra- 
tio for both CS and bremsstrahlung radiation as a function 
of the spectral index of the CS radiation, a, for m — 10 
and 10*. To enable direct comparison of those two com- 
ponents, we show the luminosity from the inner region 
(with constant B) only, assuming ro = 20 and 5 = 0.2. 
Thick and thin curves give Lcs and Lcb [see equation 
( p^ with r — ro], respectively. The relative importance 
of bremsstrahlung increases with increasing mass and size, 

e.g. Lcb/LcS OC ^0.46.-0.09 ^0.91.+0.82 ^ ^0.37^1.73 fo^ 

a ^ 0.9 and ^ 1.1, respectively. Note that Lcb decreases 
with increasing at a constant a, which is an effect of a 
strong dependence of Lcb on tt (increasing with decreasing 
G) and a weak dependence of Lcb on temperature. Also, an 
important conclusion from Figure ^ is that Eddington ratios 
^ 0.01 can be obtained from this process only for very hard 
spectra and high electron temperatures, with this constraint 
being significantly stronger in the case of AGNs. 

Figure ^ compares the upper limits on Lcs/Le with 
values of L/Le inferred from observations for a number 
of objects and with the corresponding luminosity from 
bremsstrahlung. The shown values of Lcs include the con- 
tributions from both the inner and the outer region of the 
flow [equation (p^)] for ro = 20 and S — 0.2. On the other 
hand, the shown values of Lcb are for emission from within 
ro only, and the expected actual Lcb will be by a factor 
of ~ 2 larger. For Cyg X-1 in the hard state, we used the 
brightest spectrum of Gierliriski et al. (1997), with a = 0.6, 
O = 0.2, and L ~ 3 x 10'^^ erg s~^ (assuming the distance 
of L) = 2 kpc, Massey, Johnson & DeGioia-Eastwood 1995; 
Malysheva 1997) and m = 10. For GX 339-4 in the hard 
state, we used a = 0.75, 9 = 0.1, L fa 3 x 10^'^ erg s"^ 
D = 4 kpc and m = 3 (Zdziarski et al. 1998). For NGC 
4151, we used the brightest spectra observed by Ginga and 
CG7?0/OSSE (Zdziarski, Johnson & Magdziarz 1996), for 
which a = 0.85, 6 = 0.1, and L « 9 x lO"*^ erg s"\ as- 
suming D = 16.5 Mpc (corresponding to Ho = 75 km s^^ 
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Figure 7. The Eddington ratio as a function of the spectral in- 
dex of CS radiation and for 3 values of electron temperature for 
masses (a) m = 10, and (b) m = 10*. The heavy and thin curves 
show the luminosity in the CS radiation and in Comptonized 
bremsstrahlung, respectively. See text for other assumptions. The 
curves are shown only for values of a corresponding to tt ^ 3, 
as for larger rx our formulae for a (Section Bl) break down. 



^ from an extrapolation of the 
3.1 X 10"" erg s"^ with 



Mpc-i), and m = 4x lO'' (Clavel et al. 1987). For NGC 4258 
we adopted L = 2 x lO''^ erg s 
2-10 keV luminosity of L2-iokGV 
a = 0.78 (Makishima et al. 1994) up to 200 keV, D = 6.4 
Mpc and m = 3.6 x lO'' (Miyoshi et al. 1995). Finally, for 
Sagitarius A* we assume an upper limit of L < lO^'^ erg s~^ 
(Narayan et al. 1998), and m = 2.5 x lO'' (Eckart & Genzel 
1996). Finally, we show Lcs for the average parameters of 
Seyfert-1 spectra, a = 0.9, 6 = 0.2 (e.g. Zdziarski 1999). 

We see that only in the case of Cyg X-1, with its rela- 
tively hard spectrum, the CS emission can contribute sub- 
stantially, at ^ 30 per cent, to the total luminosity. Given 
uncertainties of our model, e.g., in the value of 5, it is in prin- 
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Figure 8. Comparison of Eddington ratios inferred from obser- 
vations (asterisks, except for Sgr A*, where an upper limit is 
shown) with the model upper limits from the CS emission (heavy 
curves labeled by values of Q and a) for a range of black-hole 
masses. Vertical arrows connect observed values with the cor- 
responding CS model upper limits, except for NGC 4258 and 
Sgr A*, for which no sufficient spectral data exist. On the other 
hand, the model curves for = 0.2, a = 0.9, corresponding to 
the average spectrum of Seyfert Is, have no corresponding ob- 
servational point. Thin horizontal lines show the Comptonized 
bremsstrahlung emission. The magnetic field energy density is in 
equipartition with (a) ions in a hot accretion flow, and (b) with 
radiation energy density times I0c/v\ in a patchy corona, see 



Sections 4.1 



and 4.2.2 



respectively. 



ciple possible that the CS process can account for most of the 
emission of Cyg X-1. On the other hand, we have adopted 
here assumptions maximizing synchrotron emission and we 
consider it more Ukely that the CS process is negUgible in 
Cyg X-1. This appears to be supported by the similarity 
between X-ray spectra and the patterns of time variability 



between Cyg X-1 and the other black-hole binary consid- 
ered by us, GX 339-4 (Miyamoto et al. 1992; Zdziarski et 
al. 1998). In the latter object, Lcs clearly provides a negli- 
gible contribution to L, as shown by Zdziarski et al. (1998). 
Furthermore, a remarkable similarity exists between X-ray 
spectra of black-hole binaries and Seyfert Is (e.g. Zdziarski 
et al. 1999). For the latter, Lcs/Le ~ 10"'^-10"'' are found, 
whereas the typical Eddington ratios of those objects are 
likely to be ~ 0.01 (e.g., Peterson 1997), i.e., 3-4 orders of 
magnitude more. Concluding, CS emission in the assumed 
hot-disc geometry is unlikely to be responsible for the ob- 
served X-ray spectra of luminous black-hole sources. To ex- 
plain those spectra, an additional source of soft seed photons 
is required. Such seed photons are naturally provided by 
blackbody emission of some cold medium, e.g an optically- 
thick accretion disc or cold blobs, co-existing with the hot 
flow. 

On the other hand, the CS process can clearly be impor- 
tant in weaker sources, e.g., NGC 4258. We have compared 
predictions of our model with the 2-10 keV spectral index 
and luminosity (see above) of this object. We have obtained 
a good fit to the data, assuming S = 0.2, Q = 0.4 and taking 
into account emission within ro = 17, as shown in Figure H. 



This yields Lcs = 1.8 x 



10*^ erg s~ 



and LcB = 4 x 10 



erg s~ . A similar spectrum from an advection disc model 
was obtained by Lasota et al. (1996). 

For O « 0.2 (typical for luminous sources), we could re- 
produce the I/2-iokcV of NGC 4258, but the 2-10 keV spec- 
trum was dominated by bremsstrahlung, i.e., much harder 
than the observed one. The relatively high temperature, 
~ 0.4, required in our model, is, in fact, consistent with 
predictions of hot accretion flow models (e.g. NY95). 

As stated above, Comptonized bremsstrahlung, with 
Lcb/Le independent of M, provides the minimum possi- 
ble luminosity of a plasma with given tt, O and r. This is 
shown by horizontal lines in Figure ^. Thus, sources with 
lower luminosities cannot have plasma parameters charac- 
teristic of luminous sources, which is indeed consistent with 
hot accretion flow models, in which tt quickly decreases with 
decreasing L (e.g. NY95). On the other hand, CS emission 
can still be important in a certain energy range above the 
turnover energy even if Lcb > Lcs- This may be the case, 
e.g., in Sgr A* (Narayan et al. 1998). 



4.2 Active coronae 

4-2.1 Equipartition with radiation energy density 

Di Matteo, Celotti & Fabian (1997b, hereafter DCF97) have 
considered a model in which energy is released via mag- 
netic field reconnection in localized active regions forming a 
patchy corona above an optically-thick accretion disc. They 
assume that the energy density of the magnetic field is in 
equipartition with that of the local radiation. 



8n 



9 



167T N{rbRg 



(37) 



where rhRg is the radius of an active blob and N is the 
average number of blobs. The right-hand-side of this equa- 
tion corresponds to the average photon density in an op- 
tically thin, uniformly-radiating, sphere. This yields values 
of B somewhat larger than those of DCF97, whose used a 
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numerical factor of l/47t in their expression for B. Conse- 
quently, our estimates of Lcs will be slightly higher than 
those obtained using the formalism of DCF97. In the case 
of an optically and geometrically thin reconnection region, 
the numerical factor above should be 9/127t. Equation ds?" 
can be rewritten as 



B 



(38) 



rb \NmJ 

Following DCF97, we will assume hereafter A'' = 10. The 
characteristic size of a reconnection region has been esti- 
mated by Galeev et al. (1979) as 

I. -1/3 



(39) 



where hdRg is the scale height of an underlying optically- 
thick disc, and Qv is the standard viscosity parameter. This 
typically yields a range of 0.01 ^ rh ^ 5, depending on the 
disc parameters and increasing with accretion rate (see,e.g. 
Svensson & Zdziarski 1994). On the other hand, the specific 
corona model of Galeev et al. (1979) can power only very 
weak coronae (Beloborodov 1999) and other mechanisms can 
be responsible for formation of luminous coronae. Thus, we 
simply assume rb = 2 in numerical examples below (which 
value was also used by DCF97). 

With the above value of B, we can derive the ratio of 
the CS luminosity from the corona to the total coronal lu- 
minosity. 



Lcs _ 650y>(77/7Vmrg) °'"«°-°-°99^-"«-°-"" 
L 



]^02-24q qO. 13-1-0. 06q(]^ _ q,-) 



(40) 



for 0.4 ^ a ^ 0.9. The analogous formula for a J; 1.1 is 



^cs 
L 



1.25 1 



(41) 



y Nmr- 

From equations (^)-(^l|), we find that, in most cases, 
the coronal model yields much lower valu es of Lcs than the 



model of a hot accretion fiow of Section 4.1, For the hard 
state of Cyg X-1 with = 0.02 and m = 10 (Section 
andforiV = 10, rb = 2, we find Lcs/L~ 1.7x 10"^, which is 
~ 20 times less than Lcs / L obtained in the hot flow model. 
Only an extremely low and inconsistent with equation ( |39| ) 
value of rb ~ 10~* would lead to Lcs = L in this case. For 
objects with softer spectra, even lower values of Lcs/L are 
obtained. For objects with low luminosities, Comptonized 
bremsstrahlung would dominate, see equation (|32), which 
applies to all coronal models discussed here with r — Nr\^. 

Thus, thermal synchrotron radiation cannot be a sub- 
stantial source of seed photons for Comptonization in active 
coronal regions with equipartition between magnetic field 
and radiation, even for weak sources. This finding can be 
explained by considering consequences of this equipartition 
in the presence of strong synchrotron self-absorption. With 
neither self-absorption nor additional sources of seed pho- 
tons, equipartition between the energy densities in the field 
and in photons leads to the net Comptonization luminos- 
ity, Lc, roughly equal to the synchrotron luminosity, Ls, 
which then implies a ~ 1. (An additional assumption here 
is the validity of the Thomson limit, e.g. Rybicki & Lightman 
1979.) However, strong self-absorption in a thermal plasma 
with parameters relevant to compact sources dramatically 
reduces Ls, i.e., Ls <^ Lc- Then, a very hard Comptoniza- 
tion spectrum, with q ~ 0, is required to upscatter the few 



synchrotron photons into a spectrum with luminosity of Lc. 
On the other hand, if a ~ 1, as typical for compact cosmic 
sources, CS photons would give only a tiny contribution to 
the actual coronal luminosity, L, and other processes, e.g. 
Comptonization of blackbody photons from the underlying 
disc would have to dominate. 

This conclusion differs from findings in DCF97 as well as 
in Di Matteo, Celotti & Fabian (1999) that this process is of- 
ten important under conditions typical to compact objects, 
in particular in GX 339-4. This discrepancy stems mostly 
from their assumption of 77 = 1 (i.e., L = Le) for calculating 
the value of B in those papers (T. Di Matteo, private com- 
munication). In the specific cases they consider, Lcs ^ Le, 
i.e., only a small fraction of the dissipated power is radiated 
away, and they postulate that the remaining power is stored 
in magnetic fields (see Section 2.2 in DCF97). However, they 
do not consider the final fate of the power supplied to the 
magnetic fields. In fact, the supplied power has to be eventu- 
ally either radiated away or transported away from the disc. 
In the former case, we would recover our result of Lcs ^ L. 
In the latter, the power in magnetic fields could be either 
converted to kinetic power of a strong outfiow or advected 
to the black hole. Those possibilities require studies that are 
beyond the scope of this paper. We only note that a physical 
realization of a strong but effectively non-radiating outflow 
in the vicinity of a luminous binary appears difficult. Sim- 
ilarly, the advection time scale of an optically-thick disc is 
probably much longer than the time scale for dissipation of 
coronal magnetic fields. Finally, we note that 7] — 1 would 
often require accretion rates much higher than those esti- 
mated from the physical parameters of X-ray binaries (see, 
e.g., Zdziarski et al. 1998 for the case of GX 339-4). 

In addition, the values of O assumed in DCF97 and 
Di Matteo et al. (1999) for luminous states of black-hole 
sources are significantly higher than those used in the exam- 
ples shown in this work (which were derived from the best 
currently-available data). Due to the strong dependence of 
Lcs on Q [see eqs. (^o|)-(^)], this also results in higher val- 
ues of Lcs. 

We note that equation ( ^o[ ) can yield Lcs ~ L for 
1. This is partly explained by most of energy density 
in Comptonized photons being then in the Klein-Nishina 
limit. Then, the energy density of photons in the Thomson 
limit is much less than the total energy density, and thus 
more comparable with energy density of self-absorbed syn- 
chrotron photons. Thus, the CS process is expected to play 
an important role in sources with O ^ 1. Note, however, that 
equations (po[)-(El|) break then down as the approximations 
of equations (|l8|) and ( |2^ ) become invalid and numerical 
calculations should be employed in that case. 

4-2.2 Dissipation of magnetic field 

As discussed by, e.g., Di Matteo, Blackman & Fabian 
(1997a), the strength of the coronal magnetic field is ex- 
pected to be higher than the equipartition value of equation 
(^tI). If an active region is powered by dissipation of mag- 
netic field and the dissipation velocity is gvA, where g ^ 1 
and «A = B / {Anriamp^e / ^■lY^'^ is the Alfven speed, the en- 
ergy density stored in the field will be c/gvA times that 
in radiation (which escapes the source with the velocity c), 
and hh/vh ~ g, where hhRg is the scale height of the ac- 
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tive region (Longair 1992). Then, for a given g = hh/ri, and 
accounting for the dependence of the radiation density on 
geometry, we have. 



9 



Stt c 4(3 + ff)7t7V(rb-Rg)2c' 
Using the approximation of equation 
9.8 X 10* r;^/^ 



(42) 



we have then 



B 



5/6 
'■b " 



G. 



(43) 



= ^1/20,5/2405/24 [^(3 g)N]^/3 

The luminosity in the CS emission is then approxl 
mately given by 

Lcs _ 312 (1 + g) 

L ~ 102-25<^ [g{3 + c,)]0.61+0.3a ^ 
q3. 4-0. 3a^0. 48-0. 76a ^jY/r;)''-39 



j^0.46a-0.09Q,0.5+0.25a^]^ , 

for 0.4 Q 0.9, and 



^cs 
L 



0.59 



(l + ff)e^'°^(7V/77)°-°g 
[3(3 + ff)]''-9im0'3V0-28a0.76' 



(44) 



(45) 



for a <; 1.1, where we accounted for the dependence of the 
source area on g. 

We note that Di Matteo et al. (1997a) adopted g = 
hh/rh = 0.1, whereas Di Matteo (1998) assumed g — 1 but 
hh/rh = 0.1. We set g = hh/vh = 0.1 in our examples below. 

The values of Lcs in this model are much higher than 
those in the previous case, but still relatively low. For Cyg 
X-1 for the same parameters as above {rj — 0.02, — 10, 
Vh = 2), we find Lcs/L ~ 1, i.e. this model can, in princi- 
ple, account for the luminosity of this object. However, for 
objects with softer spectra, we obtain Lcs/L <C 1, as shown 
in Figure ib (which shows comparison with the same data 
as Figure Sk, and r; = 0.01 is assumed for the model curve 
corresponding to the average Seyfert-1 spectrum). We see 
that the predictions of this model for luminous sources, with 
r] 0.01, are qualitatively similar to those of the hot-fiow 
model. 

On the other hand, this model predicts values of Lcs for 
weak objects much lower than those in the hot disc model. 
This can be seen by comparing the dependence on rj in both 
cases. In the case of a hot flow, Lcs does not depend on 
L, so formally Lcs/L cc rj^^ (not including dependences 
of the plasma parameters on L). On the other hand, equa- 
tions (p^-(^) yield roughly Lcs/L oc r;""'^ in the present 
case (which difference is mainly caused by the source area 
much smaller in the patchy corona model than in the hot 
flow model). Thus, decreasing L leads to only a marginal 
increase of Lcs/L. In contrast, the luminosity from Comp- 
tonized bremsstrahlung is about the same in the two models. 
Consequently, Lcb 3> Lcs is typical of weak sources in the 
coronal models. 

This is indeed the case for NGC 4258, where we have 
found it impossible to reproduce its 2-10 keV power-law 
index and luminosity even with rb <C 1, except for B > 
1. Although such a high temperature cannot be ruled out 
observationally at present, the corresponding tt ^ 1, which 
then would lead to a curvy spectrum reflecting individual 
scattering profiles. Then, the similarity of the X-ray spectral 
index of that object to the average index of Seyfert Is would 
have to be accidental, which we consider unlikely. 



0.1 



0.01 



0.001 



0.0001 



10 





1 1 1 1 1 1 1 1 1 1 1 1 _ 


\ 

\ 


rj-lO-s 




'q-]0-' 


\ ^ \ 








\ ^ 




\ ^^X^ 

\ \ ' 
\ \ 

\ \ \ 

" \ \ 
\ \ 

\ '« 
\ 


\ — 










- L,g, m = 

r hs' ^ = 

= ,,,!,,, 


10» ^^^^^^^ : 

10 ^^^^C^ 

, , , 1 , , , 1 , , , : 



0.4 0,6 0.8 1 1.2 1.4 



Figure 9. The fraction of the total coronal luminosity produced 
by the CS process (heavy and medium curves for m = 10 and 
m = 10*, respectively) and by Comptonized bremsstrahlung (thin 
curves applying to both values of m) as a function of the spectral 
index of the CS radiation, for r] = 10"^ (solid curves) and rj = 
10"* (dashed cu rves') in the case of magnetic field of equation 
(|42[). See Section 4.2.^ for other parameters. 



Figure ^ shows a comparison of CS and Comptonized 
bremsstrahlung emission as functions of a for two values of 
m and two values of -q, and for g = 0.1, A'^ — 10, = 2, 
O — 0.2. As expected, bremsstrahlung dominates at low 
values of rj and high values of a, and its Eddington ratio is 
independent of m. 



4-2.3 Magnetic field in the disc 

The magnetic field strength in a corona is limited by equa- 
tion (^2|) provided most of the energy of the field is dissipated 
there, which seems to be a likely assumption. However, that 
field depends on the factor g, whose value appears uncer- 
tain, and which could, in principle, be <g 1. Furthermore, 
we cannot rule out the presence of a stronger, more per- 
manent field, which would dissipate only partially. In any 
case, the coronal magnetic field strength will be lower than 
that inside the disc, Bd- The latter is in turn limited by 
equipartition with the disc pressure, as upward buoyancy 
forces will rapidly remove any excess magnetic fiux tubes 
from the disc (Galeev et al. 1979). Such an equipartition 
field was adopted by Di Matteo (1998) in calculations of 
the rate at which energy was supplied to a corona. On the 
other hand, if the viscosity is provided alone by the mag- 
netic field, its pressure was estimated as Qv times the disc 
pressure by Shakura & Sunyaev (1973), which would ap- 
proximately correspond to multiplying equations (^)-(^^ 
below by al^'^ . Numerical simulations of Hawley, Gammie & 
Balbus (1995; see Balbus & Hawley 1998 for a review) also 
confirm that the magnetic field in the disc always remains 
below the equipartition value. 

Here, we utilize formulae for the vertically-averaged disc 
structure of Svensson & Zdziarski (1994), which take into 
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account that a fraction, / < 1, of the energy dissipated in 
the disc is transported to the corona (see a discussion on 
the transport mechanisms in Beloborodov 1999). However, 
we recalculate the disc structure taking into account the 
additional pressure provided by the equipartition magnetic 
field. For consistency with the rest of this work, we assume 
equipartition with energy density rather than with pressure, 
which assumption has only a minor effect on our results. 
Expressions for below are computed at r = 11, which 
is approximately the radius at which a radiation-dominated 
region first appears with increasing accretion rate (Svensson 
& Zdziarski 1994). For lower radii, somewhat higher values 
of Bd are obtained, but the contribution of those radii to the 
total dissipated power is small. In a disc region dominated 
by the gas pressure, we have then. 



1.1 X 10°m 



(Qvm)9/20(l-/)l/2 



(46) 



where m is the dimensionless accretion rate, m = Mc? /L-e- 
In a disc region dominated by radiation pressure, we have, 

Bd « 5.2 X 10^ [Qvm(l - f)r^''^ G. (47) 

Note that the radiation-pressure dominated region appears 
only for accretion rates higher than certain value merit . Here- 
after, we will assume Ov =0.1 and that half of the accreted 
energy is dissipated in the corona, i.e., / = 1/2. We as- 
sume the efficiency of cooling-dominated accretion in the 
Schwarzschild metric, which gives the power dissipated in 
the corona of L = 0.057m/LE. 

We first compare the maximum strength of coronal 
magnetic field of equation ( ^ ) with the disc field of equa- 
tions (^-(0). For parameters relevant to compact objects, 
the disc field is always larger than the maximum coronal 
magnetic field, up to two orders of magnitude, except for 
extremely low accretion rates, m ^ 10"^''. This shows the 
self-consistency of the coronal models of Sections 

As discussed above, given the uncertainties about the 
mechanism of magnetic field reconnection in active regions 
above accretion discs, we cannot rule out the presence of 
an average coronal field strength, (_B), str ong er than that of 
equation (|42|), and limited by equations (^6^- (^^ , 



{B) = eBd, 



(48) 



where e < 1 accounts for an inevitable decay of the field 
during a fiare event. The dissipation rate of the magnetic 
field in a reconnection event is (x B^ (see e.g. section 12.4 
in Longair 1992) so the field decays exponentially. Since the 
above average is weighted by the luminosity of an active 
region, which quickly decreases with decreasing B, e can be, 
in principle, of the order of unity. Thus, we adopt e = 0.5 in 
numerical examples below. 

Assuming the magnetic field strength given by equa- 
tion (^), we obtain for disc regions dominated by the gas 
pressure, 



Lcs _ 200 X IQ-^- 

~ Q,0. 13+0. 06a (•J^ 



,0.82+0.41(1 



_ J^0.09+0.05aj^0.41a-0.18 



(49) 



0.1 



0.01 



- m-3. 0.1. 



0.001 



0.0001 



-m 10, 0.2. 



- x_. 

1=10^ = 0.2, a=U.y \ 

\ 



V . U.K 339-4 — 



\ 



\ 



m = 4xl0'', = 0.1, Q=0.^ \ 

i \ 

,1 ■ ■ 1 1.11,1 1 I I iiml > , . liXl , \ iMiil 
0.001 0.01 0.1 1 10 



Figure 10. The ratio of the CS and bremsstrahlung luminosities 
(heavy and thin lines, respectively; the latter ratio is indepen- 
dent of rh) to the power dissipated in a patchy corona as func- 
tions of the accretion rate in the case of the coronal magnetic 
field strength equal to 0.5 of that inside the disc. The solutions 
corresponding to dominance of the gas pressure and radiation 
pressure are shown for rh below and above ~ merit > respectively. 
The ch osen v alues of m, ©, a are the same as in Figure ^. See 
Section 4.2.3 for details. The vertical bars mark our estimates of 



the accretion rates in 3 cases. 



Lcs 



: 0.052 



g2.73jv^2Q3.66^0.09 



(50) 



for 0.4 ^ Q ^ 0.9 and 



for a ^ 1.1. Similar expressions can be derived for regions 
dominated by the radiation pressure. 

Figure shows the above ratios for the same values 
of a, B and m (except that m = 10* was assumed for the 
average Seyfert 1 spectrum) as in Figure ^ as functions of m 
for e = 0.5, rb = 2, and N = 10. We see that the obtained 
luminosity is comparable to the one dissipated in the corona 
only in the case of Cyg X-1. In other cases, the CS process 
gives a small contribution, comparable to or lower than the 
contribution of bremsstrahlung. 

In the case of NGC 4258, we can reproduce the 2-10 
keV power-law index and luminosity with O — 0.43 and Tb — 
2. We note, however, that the CS luminosity is a sensitive 
function of the size of active regions in this model, Lcs oc 
Tb (which dependence is much stronger than that in the 
previous coronal models). Based on equation (^9|), we expect 
rb decreasing with the decreasing L. If this is indeed the case 
for NGC 4258, the model of coronal CS emission could be 
then ruled out for t his ob ject (unless O > 1, as in the case 
discussed in Section 4.2.2 above). 

We therefore again conclude that, in coronal models, the 
CS radiation can be important only in the case of stellar- 
mass sources with hard spectra. This radiative process is 
negligible in the case of objects with soft spectra (probably 
including low-luminosity sources) when Comptonization of 
photons from a cold accretion fiow is expected to dominate. 
Note that if the disc magnetic field pressure is less than that 
of equipartition, e.g., by (Shakura & Sunyaev 1973), the 
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predicted upper limits of Lcs would be even lower (e.g. by 
an order of magnitude for Qv — 0.1). 

We can compare our conclusions with those of Ghis- 
ellini, Haardt & Svensson (1998), who also consider mag- 
netically dominated, patchy coronae and compare the rela- 
tive importance of Comptonization of synchrotron and disc 
photons. They, however, do not assume a thermal electron 
distribution but instead assume a given form of electron ac- 
celeration and then calculate the electron distribution tak- 
ing into account a balance between the acceleration, escape, 
Compton cooling and synchrotron emission and reabsorp- 
tion. A resulting electron distribution consists typically of 
a quasi-thermal hump and a high-energy tail, which is an 
important difference with respect to the calculations pre- 
sented here. Furthermore, they assume a constant magnetic 
field for a series of models with varying luminosity, and their 
magnetic field is not directly compared to that produced in 
the disc. 

Still, their findings are compatible with ours. Namely, 
among the specific models they present, the CS process dom- 
inates only when O ~ 1 and the magnetic energy density is 
more than one order of magnitude higher than that of ra- 
diation, see their Figure 4. We find that the parameters of 
their models with the CS process dominating, correspond 
approximately to ~ 10~*, i.e. to low-luminosity sources, 
in our model above (assuming rb ~ 1 and N ~ 10). Then, 
at ~ 1, we obtain predictions about the significance of the 
CS process similar to theirs. 



5 APPLICATIONS TO ACCRETING 
NEUTRON STARS 

Power-law X-ray spectra are, in fact, observed not only from 
black-hole sources. Similar power laws are detected, e.g., 
from binary systems with weakly-magnetized neutron stars 
in the so-called low spectral state. Unambiguous classifica- 
tion of some of those sources as neutron stars comes from 
detections of type-1 X-ray bursts from them. Those spectra 
are often modelled by thermal Comptonization (e.g. Barret 
et al. 1999, hereafter B99), in which case the question arises 
of the origin of seed photons. 

It is also of interest that accreting weakly-magnetized 
neutron stars show a number of similarities to accreting 
black holes (see e.g. a discussion in B99). First, this con- 
cerns timing properties in X-rays, namely time lags between 
soft and hard X-rays, which suggests similar geometries of 
the sources and/or similar variability mechanisms (e.g. Ford 
et al. 1999). Second, spectra of black holes in their hard 
states fit remarkably well a correlation between the 1-20 
and 20-200 keV luminosities seen in neutron-star binaries 
(Barret, McClintock & Grindlay 1996) which may indicate 
that similar radiation mechanisms are operating in the two 
classes of objects. On the other hand, the presence of the 
neutron star results in more sources of soft photons than in 
the case of a black hole, namely thermal radiation from the 
surface of the star and synchrotron radiation in the stellar 
magnetic field. 

Cyclotron radiation as a source of soft photons Comp- 
tonized in weakly magnetized accreting neutron stars has 
recently been considered by Psaltis (1998). In his model, self- 
absorbed cyclotron photons are produced in a layer above 



the stellar surface with kT ~ 20 keV and a magnetic field of 
B ^ 10^" G, and then subsequently Comptonized in a hot 
spherical corona of kT ~ 10 keV. Here, we consider a more 
generic geometry in which synchrotron photons are pro- 
duced and Comptonized in the same medium, the tempera- 
ture of which is determined from observations. This leads us 
to investigating plasmas of temperatures higher than those 
considered by Psaltis (1998). As in Section 4, we compare 
the luminosities in our model spectra to those observed. 

We analyze here spectra of four X-ray bursters. The 
first one is 4U 0614-f091, which low-state spectrum was fit- 
ted with the Comptonization model of Poutanen & Svensson 
(1996) by Piraino et al. (1999). They found the seed black- 
body photons have kTsB ^ 30 eV, and the plasma param- 
eters are given by ~ 0.5 and a = 1.44. The luminosity 
was I/i-200kcV — 3.7 x 10^^ erg s~^ assuming D = 3 kpc 
(Brandt et al. 1992). It is interesting that the above elec- 
tron temperature is much higher than that seen from other 
X-ray bursters so far. The parameters of the second object in 
our sample, IE 1724-3045, modeled by B99 with a Comp- 
tonization model of Titarchuk (1994) are /cTbb — 1 keV, 
e ~ 0.057, a = 0.95. (Since B99 did not give the spectral 
index of their fit, we have obtained it from the Comptoniza- 
tion model they used.) The luminosity, given D — 6.6 kpc 
(Barbuy, Bica & Ortolani 1998), is Li_2ookoV ^ 1.3 x lO^'^ 
erg s"\ The third source is GS 1826-238, for which B99 
found a = 0.73, and Q ~ 0.08 from the high-energy tail 
with the model of Poutanen & Svensson (1996). At -D = 7 
kpc (B99), Li_200keV — 1-5 x 10''® erg s"'^. Observations of 
the fourth object, XB 1916-053, were reported by Church 
et al. (1998). The X-ray spectrum was fitted as a power-law 
of index a = 0.61 and the temperature was estimated from 
the cut-off in the spectrum to be O ~ 0.06. At D = 9 kpc, 
l'0.5-200kcv ^ 1-1 X 10^*^ erg s"^ (Church et al. 1998). 

We consider then a generic model with a spherical, uni- 
form, plasma cloud, and investigate the relation between the 
size of the cloud and the magnetic field strength required to 
reproduce the observed luminosities at their spectral indices 
and temperatures. From equations (^ and (|ll), we expect 
that this relation will be roughly of the form of B cx R~°'^ . 
We note that since a CS spectrum is normalized by its self- 
absorbed, optically-thick, part, the luminosity is roughly 
proportional to the source area, with the source geometry 
being only of secondary importance. Thus, our results may 
also approximate those for emission from a thin layer above 
the stellar surface. A second constraint on the model can 
be obtained from the fitted temperature of the seed pho- 
tons, which is known for 4U 0614-1-091 and IE 1724-3045. 
In those cases, we can identify the seed-photon temperature 
with the maximum possible turnover energy, hut ^ kTsB, 
which then yields an upper limit on B. This limit is B ^ 10^ 
G and B ^ 6 X lO** G, respectively. In the case of GS 1826- 
238, the seed-photon temperature was not fitted, but we can 
estimate, from the shape of the spectrum presented by B99, 
that fcTBB ^ 1 keV, which yields B ^ 5 x 10^ G. We have 
no corresponding limit for XB 1916-053. 

The size-magnetic field relation obtained for the four 
objects is shown in Figure |l^. Those results suggest that 
Comptonization may take place in a corona of the size com- 
parable to the stellar radius (~ 10® cm) . The required mag- 
netic field is then B ^ 10^ G. This value is consistent both 
with the upper limits resulting from fitting the energy of 
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Figure 11. The relation between the size of the emission region 
and its magnetic field strength required for the CS process to 
reproduce the observed X-ray spectra and luminosities of four 
X-ray bursters (heavy lines). The thin horizontal lines show the 
upper limits on B resulting from fitted constraints on the energy 
of seed photons. Both constraints have to be satisfied for the CS 
process to dominate. 



seed photons and with the maximum field strength possible 
in X-ray bursters, B ^ 10^°-10^^ G, see Lewin, van Paradijs 
& Taam (1995). (Their constraint results from the require- 
ment that the magnetic field does not funnel matter towards 
the magnetic poles, and more stringent constraints can be 
obtained for a specific structure of the accretion flow.) 

We note that if Comptonization takes place close to 
the stellar surface, reflection albedo close to unity would be 
required. Otherwise reprocessing in the surface layers would 
result in a strong blackbody component, which is not seen. 

We also find that it is unlikely that the size of the region 
where synchrotron photons (dominating the supply of seed 
photons) are produced and Comptonized is much larger than 
the neutron star radius. Such a configuration would be then 
similar t o th at of an optically thin disc around a black hole 
(Section 4.1) or an advection-dominated extended corona 
above a cold disc (Narayan, Barret & McClintock 1997). 
Then the required magnetic field, B ~ 10* G at f? ~ 10^ 
cm, would be too strong to be sustained by the disc [see 
equation (jsij), where we estimate the magnetic field in hot, 
optically thin, discs]. 



6 DISCUSSION 

A general feature of models of accretion fiows around black 
holes considered above is the dependence of _B oc m"^^'^ [or 
very close to it, equation (p6[)]. This then implies decreas- 
ing values of Lcs/Le with m, and the relative importance 
of the CS process much weaker in AGNs than in black-hole 
binaries. The predicted values of Lcs strongly increase with 
increasing spectral hardness and temperature of the plasma. 
Compared with observations, we find that both the hot flow 
and coronal models can, in principle, explain emission of 
the hardest (q ~ 0.5) among luminous (Eddington ratios of 
~ 0.01), stellar-mass, sources in terms of the CS process. 
However, this process provides a negligible contribution to 
the luminosity, L, of luminous stellar-mass sources with soft 
spectra and of luminous AGNs regardless of a. Taking into 
account the overall similarity between properties of luminous 



sources with either hard or soft spectra and containing ei- 
ther stellar- mass or supermassive black holes (e.g., Zdziarski 
1999), it is likely that the CS process is in general energeti- 
cally negligible in luminous sources. 

However, the hot flow and coronal models differ in their 
predictions for the relative importance of the CS process 
with decreasing L/Le- In the hot flow model, there is clearly 
a range of L/Le <^ 0.01, in which the CS process can yield 
a dominant contribution to L even in the case of AGNs. On 
the other hand, the coronal models have rather weak depen- 
dences of Lcs/L on L/Le- However, at some low value of 
L/Le, bremsstrahlung emission becomes dominant. There- 
fore, we have found that, in the coronal models, there is no 
range of L/Le (except for hard, stellar- mass, sources, see 
above) in which the CS process could dominate energeti- 
cally. In both models, bremsstrahlung is expected to give 
the main contribution to L in very weak sources. 

There are three complications in the picture above. 
First, the plasma parameters are expected, in general, to de- 
pend on L/Le, which will modify theoretical dependences of 
Lcs / L on the Eddington ratio. This is especially the case of 
the electron temperature, which is expected to be higher in 
low-luminosity sources, due to less efficient plasma cooling. 
This would strongly increase the efficiency of the CS process. 
We did not consider temperatures much higher than those 
observed in luminous objects, thus our conclusions about the 
role of the CS process in low- luminosity sources are rather 
conservative. Unfortunately, we have as yet no observational 
constraints on the plasma temperature in such sources. 

Second, we have considered here accretion onto non- 
rotating black holes. Black-hole rotation increases, in gen- 
eral, the efficiency of accretion, and the CS process may 
yield then luminosities higher than those found here. This 
is, in fact, confirmed by a study of Kurpiewski & Jaroszynski 
(1999), who have considered advection-dominated flows and 
have found that an increase of the black-hole angular mo- 
mentum leads to increase of the CS emission. This happens 
due to both the synchrotron emission itself as well its Comp- 
tonization becoming much more effective in innermost re- 
gions of the disc. Analysis of those issues is, however, be- 
yond the scope of this work. Third, even when other pro- 
cesses dominate energetically, there can be a range of pho- 
ton energies in which the CS process gives the dominant 
contribution. 

Still, the strong conclusion of Section ^ is that the ther- 
mal CS process does not dominate observed X-ray spec- 
tra of luminous black-hole sources. Rather, the dominant 
X-ray producing process in those sources appears to be 
Comptonization of blackbody photons emitted by some 
optically-thick medium, e.g., an optically-thick accretion 
disc or clumps of cold matter. Observationally, the presence 
of such component is indicated by the wide-spread pres- 
ence of Compton reflection from a cold medium in lumi- 
nous black-hole sources (e.g. Zdziarski et al. 1999) as well 
as observations of blackbody-like soft X-ray components in 
spectra of black-hole binaries (e.g. Ebisawa et al. 1994). Fur- 
thermore, a very strong correlation between the reflection 
strength and the X-ray spectral index is present in those 
sources (Zdziarski et al. 1999), in particular in Cyg X-1 and 
GX 339-4 (Gilfanov, Churazov & Revnivtsev 1999; Revnivt- 
sev, Gilfanov & Churazov 1999). The presence of such corre- 
lation is naturally explained in models with plasma cooling 
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due to blackbody emission of cold matter, whereas it cannot 
be explained if the CS process dominates. 

On the other hand, we find in Section ^ that the ther- 
mal CS process can easily account for power law emission 
of weakly magnetized neutron stars in their low states. A 
possible caveat for that model is that a correlation between 
reflection and a has been observed in 4U 0614+091 (Piraino 
et al. 1999), although its interpretation remains ambiguous 
due to the lack of a detection of blackbody emission strong 
enough to provide seed photons for Comptonization. 

We point out that a small non-thermal tail in the elec- 
tron distribution due to acceleration (as in the model of 
Ghisellini et al. 1998), e.g. stochastic or in reconnection 
events could significantly increase the turnover frequency 
(since electrons with 7 S> 1 are usually responsible for emis- 
sion at the turnover frequency), and thus increase the CS lu- 
minosity from accreting black holes and neutron stars. This 
issue is the subject of our work in progress. 

Finally, we notice thermal Comptonization may be 
(at most) a secondary process in spectral formation in 
some classes of compact X-ray sources. One important class 
of such sources are black-hole binaries in the soft state, 
which X-ray and soft 7-ray spectra appear non-thermal (e.g. 
Gierliriski et al. 1999). In that case, the X-ray spectrum 
contains a very strong blackbody component, which prob- 
ably dominates over synchrotron photons as seeds for non- 
thermal Comptonization. 



7 CONCLUSIONS 

The main results of this work can be outlined as follows. 

We have derived and tested analytic approximations for 
the synchrotron emission coefficient in thermal plasmas, ap- 
plicable especially to the semi-relativistic range of temper- 
atures, T ~ 10^ K. We have also obtained analytic approx- 
imations to the turnover frequency (at which the plasma 
becomes optically-thick to absorption). 

Then, we have presented a method to treat thermal 
Comptonization of synchrotron radiation. Our approximate 
analytic expressions allow for self-consistent calculations of 
the spectrum and the luminosity of such a source and can 
be easily applied to models of accretion discs. 

We have also investigated the role of the thermal CS 
process in accretion flows around stellar-mass and super- 
massive black holes. We have considered two main scenar- 
ios: a hot, two-temperature, optically-thin flow and active 
regions above a cold, optically-thick disc. We have found 
that this process is only marginally important in luminous 
X-ray sources containing accreting black holes, and it can 
possibly dominate only in stellar-mass sources with hardest 
spectra. The dominant radiative process in those sources 
appears to be Comptonization of blackbody radiation emit- 
ted by cold matter in the vicinity of the hot plasma. On 
the other hand, the CS process can explain X-ray spectra 
of weaker sources, e.g., low-luminosity AGNs, but only in 
the hot-flow model. Finally, below certain low luminosity, 
bremsstrahlung becomes the dominant process. 

Finally, we considered the case of weakly-magnetized 
accreting neutron stars. We have found that their power-law 
X-ray spectra in the low state can be accounted for by the 



CS process taking place in a corona of the size comparable 
to the stellar radius. 
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